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Abstract 

The  problem  of  determining  the  zeros  of  the  Hankel  function  II  (x) 
for  fixed  real  x  and  variable  (conrplex)  v  plays  a  role  in  propagation  theory 
and  in  the  theory  of  diffraction  by  a  metallic  sphere.  The  equation  li  (x)  =  0 
defines  v  as  a  function  of  x.  The  qualitative  behavior  of  this  function 
(for  X  >  O)  is  determined  for  all  x  >  j-  .      Complex  values  of  x  are  also  con- 
sidered.  In  the  case  of  a  sphere  of  finite  conductivity,  the  corresponding 
transcendental  equation  is  similarly  treated. 

Table  of  Contents 

Page 

1.  Introduction  1 

2.  The  order  of  growth  of  r(z)  with  respect  to  v  2 
5.  Zeros  of  H^(z)  for  a  fixed,  complex  z  8 
h.     Roots  of  H^  ^(z)  =0  for  real,  positive  Tallies  of  z  I5 

5.  Behavior  of  the  zeros  as  x  ->  0  17 

6.  The  derivative  of  v  19 

7.  Behavior  of  the  large  zeros  S**- 

8.  The  zeros  of  the  Wronsklan  28 

9.  SpecieJ.  ceises  35 
Appendix  3^ 
References  38 


-1- 


1.   Introduction 

The  theory  of  diffraction  of  electromagnetic  vaves  by  a  sphere  requires 
a  knowledge  of  the  zeros  of  certain  transcendental  functions. 

In  the  simplest  case  of  a  perfectly  conducting  sphere  of  radius  a, 
the  equation  in  question  is  H  (ka)  =  0,  where  a     is  the  ■Hankel  function  of  the 
first  kind  of  order  v,  and  k  is  the  wave  number  of  the  incident  wave.  The 
quantity  to  be  determined  is  the  order  v,   with  ka  being  regarded  as  a  parameter. 
Special  attention  is  paid  to  the  case  in  which  ka  is  real. 

Essentially,  then,  the  question  is:   Given  z,  for  which  values  of 
V   is  E^(2)  =  0?  While  we  do  not  answer  this  question  explicitly,  we  do  investi- 
gate the  behavior  of  v  as  a  function  of  the  argument.   Specifically,  we  derive 
certain  inequalities  satisfied  by  the  real  and  imaginary  parts  of  v  and  the 
argument  z;  we  investigate  the  behavior  of  v  as  z  — )  0;  we  obtain  a  workable 
expression  for  -r-,    and  we  find  the  approximate  location  of  the  large  zeros. 

A  more  complicated  function  must  be  considered  when  the  conductivity 


of  the  sphere  is  finite.   This  function,  W  ,  is  the  Wronskian  of  a     end  J   .   The 


V 

argianent  of  BT  depends  on  the  ambient  dielectric  constant,  and  that  of  J  on  the 
dielectric  constant  of  the  sphere.  The  techniques  used  in  the  case  of  the  per- 
fectly conducting  sphere  are  not  so  applicable  here  because  much  less  is  known 
of  W  than  of  H  .  However,  we  do  find  that  the  large  zeros  of  the  Wronskian 
behave  in  the  same  manner  as  those  of  the  Hankel  function.   In  addition  we  obtain 
a  series  of  conditions  under  which  W  ,  does  not  vanish. 


Although  many  of  the  results  can  be  extended,  the  argument  z  of 


H^ 

V 

will  be  restricted  for  the  most  part  to  the  right  half  plane;  i.e.,  |arg  z|<  p  • 
This  represents  the  physically  significant  case. 
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As  usual,  a  bar  will  denote  the  conjugate  complex  of  a  quantity. 

Because  of  the  relation  H-(z)  =  H  (z),  any  result  concerning  a  zero  of  H  can 
be  interpreted  in  terms  of  a     ;   e.g.,    for  real  z,   the  zeros  of  II  and  H  are 
conjugates . 

Another  immediate  resiilt  is  the  fact  that  the  zeros  of  n     are 
symmetric  about  the  origin,  since  H  (z)  =  e    ti  (z),  so  there  is  no  loss  of 
generality  in  letting  the  real  part  of  the  zeros  be  non-negative,  as  ve  shall 
sometimes  do. 

Throughout  this  paper,  the  letters  x  and  y  will  be  reserved  for  the 
real  and  imaginary  parts  of  z,  and  a  +  i^  will  likewise  denote  v. 

2.  The  order  of  growth  of  IT"  (z)  with  respect  to  v 

The  preceding  results,  and  all  the  subsequent  ones  predicated  on  the 
vanishing  of  M  ,  would  be  true  vacuously  if  there  were  no  zeros.  To  lend  sub- 
stance to  these  results,  therefore,  it  is  necessary  to  prove  the  existence  of 
zeros.  For  this  purpose  we  shall  apply  function  theoretical  considerations. 
Specifically,  since  E  is  not  only  analytic  but  entire  in  v,  we  shall  apply  the 
concept  of  the  order  of  an  entire  function  to  prove  that  there  exists  an  infinite 
nianber  of  zeros,  and  also  to  obtain  an  infinite  product  representation. 
From  Heine's  formula  (see  flj,  p.  26),  for  0  <  arg  z  <  «, 
g(v)  =1  e(^  -^  l^«i/2  H^Cz)  =  J^e^^  ^°^^  *cosh  vt  dt. 

If  we  let  I V I  =  p, 

\    r    w   ^     (       -y  cosh  t  +  pt  , . 
ie(^)l  <  Jo  ®  "^  dt. 

We  shall  show' that  the  order  of  g(v)  is  <  1;  i.e.,  that 
1  +  € 
|g(v)|  <  Me*^      for  all  e  >  0,  as  p  — ^  oo  , 

where  M  is  a  suitable  quantity  which  may  depend  on  e  but  is  independent  of  p. 

It  suffices  to  show  that 


\-y   cosh  t  +  p(t  -  p^)  ^^  <  j^_ 
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^     ^^   -y  cosh  t  ,,   f°°  -y  cosh  t  +  t^^^^^^/^^  • 

se  integr 


But 

I 

—       in 

P 

Since  each  of  these  integrals  converges. 


e(v) 

for  large  p;  i.e.,  the  order  of  g(v)  is  <  1. 

2  1/2 

Now  let  V  =  X.  Then  g(v)  =  g(X  '  )  s  f(x)  ig  an  entire  function 

of  X.  Moreover,  since 

the  order  of  f(X)  is  <  1/2.  Hence  (see  [2],  pg.  250), 
f(x)  =  f(0)  TT  (1  -  ^), 

where  the  X^  are  the  zeros  of  f(X).  We  note  that  f(0)  /  0  since  ti  (z)  ^  0 
from  Corollary  3.3^  to  be  proved  later. 

In  terms  of  v,  this  product  hecomes 

hj  -  e-^/=  ^  u  (^  -  y 

where  the  v,  are  the  zeros  of  n  (z)  for  y  >  0.  Moreover 

f(x)  =  f(o)  IT  (1  -^)  =  f(o)  [1  -  ^^   r^  ••■] 

=  f(0)    +  X  f'(0)    +   ... 

and  we  find 

^    ^  =  -f'(0)/f(0). 


x:  3^     siiv)  dv  jo  t%^^  "°"^  ^dt 

k        2-       nVy     dX       v=0--      ^   izcosht,, 
k  2j     e  dt 


Since  f(x)  is  an  entire  function  of  order  <  1  and  is  not  a  polynomial, 
it,  together  with  the  Hankel  function,  assumes  each  value  infinitely  often.   In 
particular,  when  y  >  0,  H  (z)  has  an  infinite  number  of  zeros  (see  ^2^ ,  ch.  8). 

We  can  extend  this  result  to  the  real  axis  by  applying  the  same  technique 
to  the  function  h(v)  defined  below.  When  x  >  0,  from  Nicholson's  formula 
([l],  p.  51), 

(2.1)  h(v)  =  2^H:J;(z)  H^(z)  =  L  K  (22  sinh  t)  cosh  2vt  dt, 

where  K  is  the  modified  Hankel  function  of  order  zero.  Then,  if  we  put  |v|  =  p, 
|ti(v)|  <  (^  |K^(2z  sinh  t)|  e^P'^dt. 
Now  for  any  v  and  z  for  whi  ch  x  >  0  ( [5] ,  p .  I81 ) , 

(2.2)  y  z)  =  J^e"^  '^"^^  ^cosh%  S  dS. 
Therefore, 

00   i  ! 

i,r  /-o   •  V  4-^  I  /  (     -2z  sinh  t  cosh  S 
K  (2z  sinii  t)  <     e  dS 

I  qV         '  I  _  Jo   I  I 

/   -2x  sinh  t  cosh  S,_, 
=    e  dS. 

Since  cosh  S  >  1  +  -;r-,   we  have  for  any  a  >  Q 

f°°  -2a  cosh  S^^  ^  -2a   f°°  -aS^  .^     JT     e"^^. 
Jo  ^         dS  <  e     J^  e     dS  =  ^  — :, 


2x  sinh  t 


and  |K  (2z  sinh  t)  |  <-^  ® 


Thus; 


^x  sinh  t 
j —   .00   -2x  sinh  t  +  2pt 

ih(v)i<4-^  ^-=== — dt. 

'^  '^  J   X  sinh  t 

To  show  that  the  order  of  h(v)  is  <  1,  it  suffi-ces  to  show  that 

I-^  f°°  exp  L-2x  sinh  t  +  p(2t  -  P^)J^^  <  j^^ 
'°  yx  sinh  t 
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Nov 


pV2      -2x  sinh  t  ,00  „        •■ua../-,J.^l+~ 

•^'^    e f         exp    |_  -2x  smh  t  +  (2t)  e 


:■<  i;  7==^*^  I 


■dt. 


^x  sinh  t        p  /2     ^x  sinh  t 


(•''  •  •) 


Since  each  integral  converges,  h(v)  =  O^e      j  for  all  positive  € ,  and  the 

order  of  h(v)  is  <  1. 

1/2 
As  before,  we  have  h(v)  =  h(X  '  )  =  P(x).   Since  cosh  2vt  is  an  even 

function  of  v,  P(x)  is  an  entire  function,  and  is  of  order  <  p.  Then  there  are 

an  infinite  number  of  zeros  of  P(>^),  and  hence  of  BT  H  .   If  z(=x)  is  real 

and  positive,  then 


hJ(x)   =  h|(x)  and  H^(e^"x)   =  -^x) , 

and  with  our  previous  results  we  have 

Theorem  2.1  If  0  <  arg  z  <  n,  there  are  an  infinite  number  of 
zeros  of  n  (z) . 
Since  P(X)  is  of  order  <  — , 

p(x)  =p(o)  TT  (1  -  ^)  . 


Az)  H^(z)  =Hi(z)  H^(z)  irfl--^ 


b^  '  o''  '  ' '       2 
kV     v^ 

for  X  >  0,  where  the  v,  are  now  the  zeros  of  E^  a     .     Moreover,  since 

h(v)  oh(0)  TTI  1   '^ 
k  .     ,^ 

=  h(o)  (  1  -  v2  x:  ^  + . . 
\        ^  \ 
2 

=  h(0)  +  ^  h"  (0)  +  .  .  .  , 
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we  have 


s-    -1-    =    ^"  (0) 

V       2  "2hj0) 


2   f°°t^  K  (2fesinh  t)   dt 


k    v: 


,00 

K  (2z  sinh  t)   dt 
)o     o 


For  z  real  and  positive,   if  v^"^     is  a  zero  of  H"^,  we  have  v  ,    =  v  ,  ' 


Hence 


k    v,  k 


1112  ■'J2T2 


k 


Tiys  ^  ji)2 


and 


Theorem  2.2 


^'^     JI72 


k      V, 


.00 

t     K     (2x  sinh  t)   dt 
Jo o_ 

K     (2x  sinh  t)   dt 
Jo     o  ^ 


H^ 


To  show  that  n     has  only  sliirple  zeros,  we  reqiiire  the  following 


well-known 


Lemma  2.1  Let  F(t)  be  positive,  continuoiis,  and  decreasing  for 

00 
t  >  0.  Then  J  F(t)  sin  t  dt  >  0. 


Proof: 


.00  00 

I  F(t)  sin  t  dt  =  51 


00 
n=o 


r 


2nn    ^(*)  ^^"  *  ^*  M  (2n.l)n^(*)  ^^"  * 


r: 


dt 


By  the  first  mean  value  theorem, 

/(2n+l)rt  /(2n+2)it 

J2nrt    ^^*^  ^^"^  t  dt  +  J(2n+i)/(*)  ^in  t  dt  =  2  F(f)  -  2  F(t"), 

for  some  t',  t"  such  that  2mt  <  t'  <  (2n+l)rt  and  (2n+l)rt  <  t"  <  (2n+2)rt.  Since 
F(t)  is  strictly  decreasing,  P(t')  -  F(t'0  >  0. 
We  generalize  this  resvilt  in 

Lemma  2.2.  Let  q(0)  =  0,  q'(t)  >  0,  q"  (t)  > '0,  p(t)  >  0, 

p'(t)  <  0,  with  p'(t)  continuous,  for  t  >  0. 
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Then 


.00 

p(t)  sin  q(t)  dt  >  0. 


Then 


Since 


and 


Proof: 
Let  X  =  q(t).   Since  q'(t)  >  0,  t  =  q""'"(x)  s  t(x) 


j7p(t)  Sin  q(t)  dt  =  (^  ^Ij^  Sin  X  dx. 


F(x)  a  ^,  >  0 


>  •  ^  «        TXr,  •• 


F.(x)  a  rr  -  pq"  <  o, 

the  lemma  is  a  consequence  of  Lemma  2.1. 

Now  we  can  prove  the  simplicity  of  the  zeros.  For  if  we  let  n  (z)  =  0 
where  x  >  0,  from  a  form\ila  of  Watson  (see  [ij,  p.  30) 


t(2)c  \   V  ^  '     oi   f  J,  /_   J  r.  j.\  -2vt  ,. 

r   (z)  3; = K  (2z  sinh  t)  e     dt. 

V  ^      ov        n  Jo  o  ^ 


Since 


Ti(2)/  \    /  n       J  tr  /^   •  1-  j.\    (       -2  z  sinh  t  cosh  s, 
H^   (z)  !t  0  and  K  (2z  sinh  t)  =    e  ds, 

V  ^   '        o  Jo 


it  stiff  ices  to  show  that 
,00  ,co 


Im 


j   j   exp  (-2z  sinh  t  cosh  s  -  2vt)  ds  dt  ^  0, 


.00  .00 

j   j   exp  (-2x  sinh  t  cosh  s  -  2at)  sin  (2y  sinh  t  cosh  s 

+  2pt)  dt  ds  /  0, 

where  we  changed  the  oirder  of  integration.  As  mentioned  earlier,  we  can  asBiime 

a  >  0  with  no  loss  in  generality.  Since  for  y  >  0  we  can  prove  that  P  >  0 


(see  Corollary  5-3  and  Theorem  ^.l),  the  integrand,  as  a  function  of  t,  satisfies 

the  requirements  of  Lemma  2.2  and  is  therefore  positive  for  each  s.   Consequently 

ShJ(z) 

— r 3fc  0  and  we  have 

ov 

Theorem  2.3.  The  zeros  of  r(z)  are  simple  when  x  >  0  and 
'y  >  0. 

3.   Zeros  of  H   (z)  for  a  fixed,  complex  z. 

In  this  section  we  shall  derive  several  inequalities  relating  the 
zeros  V  and  the  complex  argimients  z.  For  the  most  part,  z  will  be  restricted 
to  the  upper  half  plane.   If  we  include  the  positive  real  axis  this  corresponds 
to  the  actual  physical  situation,  since  the  wave  number,  which  is  essentially 
z,  lies  in  the  first  quadrant.   If  z  is  in  the  lower  half  plane,  the  zeros 
behave  entirely  differently.  This  is  not  too  surprising,  since  the  asymptotic 
expansions  or  the  integral  formulas  for  n     reveal  a  basic  difference  in  the 
two  half  planes.  For  instance,  when  y  <  0  and  x  /  0  it  has  been  shown  that  the 
total  number  of  z's  which  satisfy  r(z)  =0  for  real  v  "is  the  even  integer 
nearest  to  v  -  p,  unless  v  -  -  is  an  integer,  in  which  case  the  number  is 
V  -  -       ([5],  pp.  511-2).  But  for  y  >  0,  thei;p  are  no  real  values  of  v  for 
which  H^''"^(z)  =0,  as  is  asserted  in  Corollary  3.5-  The  following  res\alt  is 
due  to  B.  Friediaan: 

Theorem  3-1.  The  zeros  of  r(iy)  for  y  >  0  are  pure 

imaginary. 

Proof:   Consider  the  differential  equation 


with  u(y  )  =  u(oo)  =  0.  This  equation  is  of  the  Sturm-Liouville  type.   It 
has  the  solution  u(y)  =7y  H^(iy).   Since  the  eigenvalues  v  are  real,  v   is 


for  V   real.  Then  the  zeros  of  Ir'(iy)  are  piire  imaginary 
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either  real  or  pure  imaginary.  The  former  is  impossi'ble,  since  by  Heine's 

formula 

■a^f^      ■»     2i  ^-v«i/2   f   -y  cosh  t   ,   .  j,  /  . 

H(iyJ= e    '      e''o      cosh  vt  dt  ^  0 

v^   "" o  n  Jo  ' 

Ht(iy) 

We  have  something  of  a  converse  in 

Theorem  $.2.  Given  P  /  0,  there  exists  a  positive  y  such  that 

Hjp(iy)  =  0. 

Proof:  For  b  real,  we  have  (see  [l] ,  p.  5^) 

r^-t  cosh  b  ^(1)(.^)  ^^  .  _  ^^  ft^'H^"''      ■ 
}o  ip  ^  sinh  prt  sinh  b 

Choose  p  =  —  .  Then  the  integral  vanishes.  Kow  H7o(it)  is  piare  imaginexy,  since 
D  ip 

-^ip(^^^  =sfp(.it)  =Hjp(it). 
Because  of  exp  (-t  cosh  b)  >  0,  H  (it)  mvist  vanish  for  some  t  >  0.   Since  b, 
and  hence  3,  was  arbitrary,  the  theorem  is  proved. 

In  this  case  we  can  show  that  |p|  >  y  by  means  of  the  following  class- 
ical argument.  The  function  r(ze  )  m    u(t)  +  iv(t)  satisfies  the  equation 

^  f;(ze*)  +  (zV*  -  v^)  Hi(ze*)  =  0. 
dt 

In  terms  of  its  real  and  imaginary  parts,  this  becomes 

u"  +  [(x^  -  y^)  e^*  -  (a^  -  p^)]  u  -  2(xye^*  -  o^)  v  =  0 

-   v"  +  [(x^  -  y^)  e^*  -  (a^  -  p^)]  v  +  2(xye^*  -  o^)  u  =  0. 
First  we  multiply  these  by  v  and  u  respectively  and  subtract,  getting 

u"  V  -  u  v"  -  2(xye^*  -  oP)  (u^  +  v^)  =  0. 
If  we  multiply  instead  by  u  and  v  respectively  and  add,  we  obtain 

uu"  +  w"  +  [(x^  -  y^)  e^^  -   (a^  -  p^)]  (u^  +  v^)  =  0. 
Integrating  these  equations  then  gives  us    .  ..... 
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Theorem  30-      If  ir"(z  e  °)    =0  nnd  y  >  0,   then 
(3-1)      1^  (xye^*  -  a3)    (u^  +  v^)   dt  =  0 


( 
"o 

and 

"o  '  o 

V?e  now  let  z  =  v.     Then  equation  (j-l)  'becomes 


(5.2)  C  [(x^  -  y'^)  e^'"^  -  (a^  -  p^)]  (^2  ^  ^2^  ^^  ^  |°°(^,2  _^  ^,2^^^  >  q_ 


a6 


r  (e^*  .  1)  (u^  +  v^)  dt  =  0. 


If  t  >  0,  then  we  must  have  o^  =  0.   But  then  a  =  0  since  P  =  y  >  0,  and  (3-2) 


o  — 
becomes 


P^  j^   (1  -  e^*)  (u^  +  v^)  dt  >  0. 


This,  hovrever,  is  impossible  if  t  >  0.  Therefore,  t  must  be  negative. 
We  have  thus  obtained  the  following  result. 

Corollary  3.1.   If  r(av)  =  0  for  a  >  0  and  P  >  0,  then  a  <  1. 
We  have  a,s  a  special  case 

Corollary  3-2.      If  H^gCiy)    =  0  for  y  >  0,    then   |pl   >  y. 
V7e  can  exploit  equations   (5'l)   and  (3 '2)   somewhat  further.     For 
instance  let  x  >  0,    and  recall  that  y  >  0.      If  t  is  siifficiently  large, 

Pi" 

X  y  e   -  cxp  >  0.  Then  for  the  integral  in  (3'l)  "to  vemish,  since  the  above 

2t 
function  is  monotonic  in  t,  we  must  have  xye  o  -clP<0.  Similarly,  if 

X  <  0,  the  inequality  is  reversed.  Finally  if  x  =  0,  from  Theorem  3.I,  a  =  0. 

This  establishes 

Corollary  3.3.  If  H  (z)  =0  for  y  >  0,  then  sgn  (op  -  xy)  =  sgn  x. 


Moreover,  if  y  =  0,  then  H  (x)  =  h(x) .  But  these,  as  linearly 
'    *^  a       a 

independent  solutions  of  Bessel's  equation,  cannot  vanish  simultaneously. 
Together  with  Corollary  3'3^  this  proves  that  if  y  >  0,  then  n{z)  /  0,  a 
result  first  obtained  by  Macdonald  ([3];  P'  511)- 
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Again  from  equation  (j.l),  if  x  a  P  =  0,  then  x  =  a  P  =  0.   By- 
considering  each  quantity  separately,  vre  obtain  the  folloving  results.   First 
ve  let  a  =  0;  then  ve  have 

Corollary  3  A  If  H^q(z)  =  0  and  y  >  0,  then  z  =  iy. 
Secondly,  letting  p  =  0,  we  obtain  MacDonald's  result.   Finally,  if  x  =  0,  we 
find  H  (iy)  =  0  implies  v  =  a  or  v  =  ip.  But  v  =  a  is  impossible,  by  Heine's 

formula,  and  we  again  obtain  Theorem  J-l- 

2    2 
If  we  now  consider  equation  (5-2)  and  assimie  x  -  y  <  0,    the  same 

2    2  2t      2    2 
reasoning  as  before  indicates  that  (x  -  y  )e  °  -  (a  -  p  )  >  0.  This  is  the 


content  of 


Corollary  3-^.   If  H^Cz)  =0  with  x^  <  y^  and  y  >  0,  then 


2       J2.    .     2         2 
a  -  P  <  X  -  y  . 

Further  results  restricting  the  location  of  the  zeros  can  be  obtained 

from  Hankel's  asymptotic  formula.  This  expression  is  valid  not  only  for  the 

magnitude  of  r  of  z  large  compared  to  that  of  v  and  1,  when  it  attains  greatest 

accuracy,  but  also  for  small  values  of  z.   If  v;e  consider  only  the  first  term 

of  the  expansion,  Hankel's  formula  is 

(3.5)    hJ(z)  -  (fj)^V<'-  ■^'(i.K). 

We  shall  find  conditions  under  which  the  remainder  R  has  magnitude  less  than 
unity,  whence  H  (z)  /  0.  First  we  shall  use  Weber's  form  for  the  remainder. 

The  derivation  of  this  formiila  for  real  values  of  v  is  reproduced 
in  Watson  ([jlj  PP-  211-2).  This  can  be  extended  to  complex  values  by  only 
slight  changes.  Without  going  through  this  lengthy  analysis,  we  merely  state 
the  following  inequality  for  R.   If 

2  |z|  =  2r  >  a  -  I  >  0,  p  >  0,  and  x  >  0, 


then 


r2a-l 


r(v^i) 
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But  it  can  easily  be  sha/n  that 
■2a  -  1 


when 


2r  -  a  +  2 


2r   /  -  2r  -  (a  +  l)(2a  -  1) 

(a  +  l)  (a  -  ^)  <  r,  and  that 


r(a  +  |) 


r(v-|) 


<^sinhp«, 


Then 


r.       -,              2r-a+-r         -1.0 
|r|  <-^   \v^       ^  I ^ sinh^ 


?  '  2r  -  (a  +  l)(2a  -  1)    Pit 
vhich  is  less  than  unity  when 

0  <  k-r^  -   2r 


2r 

(a  +  l)(2a  - 

1  ^+  - 

v2 

1 
-  IT 

sinh  P« 

•  ±;  + 

P 

(a- 

•1) 

2       1 
V     -5: 

sinh  Pit 

This  is  a  quadratic  in  r,  giving  us  a  simpler  inequality  in 

Theorem  $A.      If  H^(z)    =  0  for  a  >  I,   P  >  0^    and  x  >  0^ 
then  2|z|   <  (a  +  l)(2a  -   l)    +   |    v^  -  ^    |  £i£|_Pi  . 

We  can  complement  this  result  insofar  as  a  is  concerned  by  considering 

another  form  for  the  remainder  due  to  Schlafli  ([jj^  P*  219).  For 
|a|  <  I  and  -  |  <  arg  -z.  <  ^  , 


the  remainder  R  of  Hankel's  formula  satisfies 


(3-M 


iRl  ^  I  ^^^^ 


2   1 
a  -^ 


cosh  Pit 


2   1 


)s  ojr  '     2r    -  |  cos  cue  |      2r 

where  c  =  1  when  y  >  0  and  c  =  |sec  arg  z|  when  y  <  0.  Then  if  H^(z)  =  0,  we 
mxist  have  |r|  =1,  which  gives  the  following  resiats. 
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Theorem  3.^3.      If  E^{z)    =  0  for    |a|   <  |  and  -  |  <  arg  z  <  2«   ^ 
then 


Or.  ^  c   cosh  grt      I      2       1 

^r  <  — r p         a     -  T- 

-      I  cos   cut|       '  4 


vhere  c  is  defined  es  above. 

Since    |a|   <  ^,   we  obtain  the  following  less  accurate  but  sinrpler 

inequality  for  y  >  0: 

.   cosh  Brt 

r  <  -1 ^^—r    . 

—   I  cos  an  I 

k.     Roots  of  H   (z)  =0  for  real,  positive  values  of  z 

Most  of  the  results  of  the  preceding  section  are  valid  only  when 
y  >  0.   In  this  section  we  shall  restrict  z  to  the  positive  real  a::d.s. 

The  first  result  will  follow  from  a  lemma  which  will  be  useful  later 


Lemma  k.l.      If  r(x)  =  0  and  x  >  0,  then 


Jx  '  v^  ' ' 


V 

t   oPJt 


Proof:  We  consider  the  self  adjoint  form  of  Bessel's  differential 
equation  for  H"^(x)  and  H_(x) .   Ck)ing  through  the  already  familiar  steps  of 
multiplying  each  equation  by  the  other  function,  subtracting,  and  integrating, 
we  obtain^ 

/: 


2       -2 

.CD     V       -    V_      jjl(^)    g2^^)     ^^    ^    ^ 


H^(t)  -rf  £{t)    -  H!(t)     ^     H^(t) 
vdt     V  V  dt       V 


r2  dt   _  l+i     prt 
I        t         rt   ^ 


i^icoB  1^  |Hj;(t) 

from  the  behavior  of  the  Hankel  ftinctions  at  infinity,  and  the  fact  that 
Since  the  integral  in  the  last  equation  converges,  we  have 


-Ik. 

Theorem  4.1.  If  r"(x)  =  0  for  x  >  0,  then  a  p  >  0. 
As  far  as  a  is  concerned,  this  tells  us  only  that  |a|  >  0.  We  can 
improve  this  inequality,  but  first  we  require  the  folloving  well-known  extension 
of  Lemma  2.1. 

Lemma  '4-. 2.  Let  J   f(t)  cos  bt  dt  exist  for  b  /  0,  where  f(t) 

is  twice  differentiable  for  t  >  0  with  f'(t)  <  0  and  f"(t)  >  0. 

-oo 
Then  J   f(t)  cos  bt  dt  >  0. 

Proof:   Integrating  by  parts,  we  obtain 

oo 
f^f(t)  cos  bt  dt  =  ^  f(t)  sin  bt    "  ^  f  ^'(*)  sin  bt  dt. 

.00 

Since    f(t)  cos  bt  dt  exists  with  f'(t)  <  0,  we  must  have  lim   f(t)  =  0. 
/O  t  -4  00 

And  since  J   f(t)  cos  bt  dt  exists  and  f(t)  is  monotonic,  f(t)  =  6{t~   )  as  t  — ^  0; 

hence  lim  f(t)  sin  bt  =  0.  Then 
t  ->  0 


J  f(t)  cos  bt  dt  =  -  ^  f  f'(t)  sin  bt  dt, 


and  an  application  of  Lemma  2.1  completes  the  proof. 
We  are  now  able  to  prove 

Theorem  4.2.  Let  Br(x)  =  0  with  x  >  0.  Then  |a|  >  x  -  -/x/2. 

Proof:  From  Nicholson's  form\ila  (2.1)  for  x  >  0, 

=  J  K  (2x  sinh  t)  cosh  2at  cos  2pt  dt. 


Re 


r  <(^)  ^(-) 


We  shall  e^ply  the  preceding  lemma  to  find  conditions  under  which  this  does  not 
vanish.  We  set 

f(t)  =  K  (2x  sinh  t)  cosh  2at   * 


=/:• 


.2x  sinh  t  cosh  8  ^^^^   2at  da 


for  t  >  0.  Then 
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^i/^\    (       -2x  sinh  t  cosh  sr_     ,      ,,    vo+ 
f'(t)  =    e  \_-2x   cosh  s  cosh  t  cosh  2at 

+  2a  sinh  2at]  ds  <  0 
when  |a|  <  x,  and 

f  (t)  =  J^e-2^  ^^^  ^  =°^^  ^  A(a,x,s,t)  ds 

where 

2     2      2 
A  =  ^x  cosh  s  cosh  t  cosh  2at  -  8clx  cosh  s  cosh  t  sinh  2at 


2 
■2x  cosh  s  sinh  t  cosh  2at  +  ka     cosh  2at. 


Let  a  >  0.   Then 


r,.  2   .,2.  .__,_2.    „_v,,  ..^  .  .  ,  i,-2"| 


A  >  M+x  cosh  s  cosh  t  -  2x{ha,+l)   cosh  s  cosh  t  +  ^a  |  cosh  2at 


which  is  non-negative  for  all  s  and  t  when 


Thiis  f(t)  satisfies  the  conditions  of  Lemma  k.2   and  we  find  that  if  H~(x)  =  0 
with  a  >  0,  then 


^  4a  +  1  +JQa,  +  1 
X  <  j^ . 

But  this  can  be  shown  to  imply  a  >  x  - Jx/2   . 

If  we  restrict  [a]  to  less  than  — ,  we  obtain  another  inequality. 

Theorem  4. 3.  Let  tt  (x)  =  0-with  x  >  0  and  |a|  <  -^  . 
Then 

X  <  P  tan  ojt  +  -^  . 

Proof:  Using  a  formula  due  to  Watson  {[}\,   p-  ^^5)^  we  find  for 

X  >  0  and  |al  <  ■^  that 

8i  Sin  2a^  r        (2^  ^.^  ^j  ^-2ipt  ^^ 
)o     2a 


2 
jt 


=  H^x)  h2(x)  -  hJCx)  Hf-(x) 

=   |H^(x)|2e^--  |Hi(x)|V^^  • 
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Then 


(if.l)   COS  cut  f  K2J2X  sinh  t)  cos  2pt  dt  =  |g  (|H^|^  e^'^  +  \E^J^   e"^'^), 
(h.o)     sin  cut  J^  K2J2X  sinh  t)  sin  23t  dt  =  ^  ([H^I^  e^  -  |hJ|^  e"^"^)  . 


After  subtracting  we  obtain 
2 


^   |H^(x)|^e'^'^  =  Ck^^  (2x  sinh  t)  cos  (2pt  +  cut)  dt. 
Dividing  both  sides  by  cos  cut  and  integrating  by  parts  give  vis 


2     -Pn|„l|2 
rt     e         H 


/r^i"2pt|  2^     ^K2j2xsinht) 


+  tan  cut  K     (2x  sinh  t)  ^  dt. 


Applying  equation  (2.2)  for  K  (z)  yields 


2  -p«  |H^|^    ,00 
jt  e  ^   '  v'     ^ 


J   f(t)8in  23tdt 

o 

where 


Then 


-/,  \    I   -2  X  sinh  t  cosh  s,^/x    ,^    ,     ^     n^ 
f{t)    =  I  e  cosh  2as  (   ^  cosh  t  cosh  s  -  tan  out)  ds. 


41.1 /■  4- ^    f       -2x  sinh  t  cosh  s   ,  ^     , 
f'(tj  =  I   e  cosh  2cus  cosh  s 


^  sinh  t  cosh  s  -  2x  cosh  t  (  —  cosh  t  cosh  s  -  tan  cut) 


ds. 


If  P  >  0  and  X  >  p  tan  cut,  then  f(t)  >  0;  also  f '(t)  <  0  when  sinh  t  <  2  cosh  t 

(x  cosh  t  cosh  s  -  P  tan  cut).  Simplifying  this  and  applying  Lemma  2.1  complete  the 

proof. 

We  can  close  the  a  interval  under  consideration  in  the  following  manner. 
If  we  let  V  =  -  +  ip,  then 
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nt  ^  V   V  -  1     V-1   V       ^     V       -V  -^       y^  ' 

where  the  argument  of  both  Hankel  functions  is  t.   This  is  easily  proved  for 

a  =  -  — ,  and  incidentally  provides  an  extension  of  C+.l).  Nov  if  we  divide  both 

sides  by  t  and  integrate  from  x  to  od  ,  applying  Lemma  ^+.1  we  find 


,00 

_       ^  J^    \ 
jtx       Prt 


-'--  '-^'"/r  i<(t)i'f ' 


whence 


Theorem  hX.     Let  Br(x)  =  0  with  x  >  0,  a  =  +  t.   . 
Then  |p|  >  |  . 

If  we  apply  Lemma  2.1  to  equation  (^-2),  we  obtain 
Theorem  ^.^.   If  x  >  0  and  |a|  <  p,  then 

sgn  [|H^(x)|  e^  -  JH^Cx)  |]  =  sgn  o^. 

5.  Behavior  of  the  zeros  as  x  — ^  0 

The  inequalities  |a|  >  x  -_/x/2  and  x  <  3  tan  cut  +  ^,  which  we  found 
in  the  preceding  section,  give  us  no  information  when  0  <  x  <  -.  We  can  partially 
explore  this  range  by  letting  x  tend  to  zero,  which  enables  us  to  employ 

approximation  formulas  of  a  relatively  simple  nature. 
Let  us  consider  Lemma  ij-.l: 

where  H'^(x)  =  0  for  x  >  0.  If  we  let  x  — ^  0,  the  integral  diverges,  whence  either 

oP  — >0  or  p  — >  00.  The  latter  possibility  is  eliminated  since  H  =0  implies 

n  e"^  '' 

n       =0;  and  if  P  — ^  00  then  g.    — ^  0  unless  a  — ^  0,  whereas  the  integral 

diverges.   Consequently  either  a  — ^  0  or  p  — )  0. 

We  can  show  that  at  least  a  ajrproaches  zero. 
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Theorem  ^.1.  Let  H  (x)  =0.   If  x  — )  0+,  then  a  — >  0. 
Proof:   Let  H  (t)  =  u(t)  +  iv(t).   Then  separating  the  self-adjoint 
form  of  Bessel's  equation  into  real  and  imaginary  parts,  we  obte,in 


and 


(tV)'   +(t.2i^)  v  =  2|g^. 


Miiltiplying  by  u  and  v  respectively  and  adding: 

2  2 

(5-1)  (tu')'    u  +  (tv')'   V  +  (t  -  "•     :   ^   )    (u^   +  v^)    =.0. 


Integrating  from  x  tc  as 

t  (u'2  +  v'^)  +  (t  - 


c 


•2  .  ,..2^  .  f.        Z^)  (,2  ,  ^2) 


dt  =  0. 


From  Lemma  ^-1,  this  simplifies  to 


P 

dH^ 

V 

dt 

Let  X  — )  0.   If  a  >  e  for  some  e  >  o,  then  as  shown  below  the  integral  diverges 

2    2 
to  -  CD  ,  and  a  <  P  .  But  since  a  — >  0  or  p  — >  0,  this  implies  that  a  — )  0. 

To  show  that  the  integral  diverges  as  stated  we  consider  the  integrand 

for  small  t,  and  a  >  e  >  0.  Then 


|Hi(t)|    = 


J     (t)    -  e"^^  J  (t) 

_  V  V 


rt    ^2' 


r(v) 


Similarly 


dt      V 


rt    ^2^ 


-a-1 


r(v  +  1) 


Thus  the  second  term  dominates  the  integrand,  and  the  integral  diverges 
to  -  00.  This  completes  the  proof. 

For  small  x,  we  now  find  a  relation  involving  both  a  and  p. 
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Theorem  ^.2.  Let  H''"(x)  =  0.   If  x  — ^  0+,  then 

lim{a  log  I  +  ^)  =  0. 

Proof: 

J  (x)  -  e-^^'^  J  (x) 


Let  X  — >  0+,  and  restrict  a  to  positive  values  also.   Then 

lim   l(|)"V(v)|  =lime-^l°Sx/2jp(^)|  ^  ^  ITCvJI 
X  -^  0+  '  I  ^ 


where 


Also 


c  >  1  and  V  =  i^  =  11m  v. 
—       o     o 


lim  |(|)V(-)e-^^|  =  lim  e-  1°S  ^/2  .  p«  \r{-.)\    =  ^  e^^''  iPC-v^)! 

Since  v  =16,  If^Cv  )|  =  IF^C-v  )l:  and  since  Er(x)  =  0  in  the  limit  also, 
o     o   I  '   ^  o  I    I  '     o  I '  V 

6  rt/2       _^  ^^  ^^/^ 

the  tn/o  limits  above  miost  be  equal,  whence  c  =  e     =  lim  e     ^     '     . 

6.  The  derivative  of  v 

We  can  obtain  some  interesting  results  not  only  for  small  x,  as  in  the 

dv 
last  section,  but  for  any  x  by  examining  -t— .  However,  it  is  not  much  more  difficult 

to  consider  complex  values  of  the  argument.  W.e  shall  derive  a  formula  which  was 

first  obtained  by  Watson  (see  [5],  p-  ^08) , 

^  =  2c  I  K  (2c  sinh  t)  e'^^^   dt, 
dv     jo  o^  .   ' 

where  v  is  real  and  c  is  a  positive  zero  (in  z)  of  the  cylindrical  function 

J  (z)  cos  s  +  Y  (z)  sin  s  for  constant  s.   This  form  immediately  excludes  H  . 

V/e  shell  generalize  it  to  apply  to  complex  values  of  z  and  v,  as  well  as  to  any 

function  of  the  form  C  (z)  =  c,J  (z)  +  c_Y  (z),  where  c,  and  c^  are  independent 
V      1  v^      2  V  1      2        ^ 

of  both  z  and  v. 

All  the  results  of  this  section  will  be  corollaries  of 

Theorem  6.1.   If  C  (z)  is  of  the  above  form,  with  x  >  0,  then 
^any  zero  v  satisfies 


(6.1)     2z  f^  J7Ko(2z  sinh  t)  e'^"*  dt  =  1. 
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Proof:   Let  B  (z)  =  b,  J  (z)  +  b^  Y  (z)  be  a  function  of  the 
V       1  V       2  V 

same  type  as  C  (z)  and  linearly  independent  of  it-   Then 
5  =  b^c^  -  bgC^  /  0. 

The  V/ronskians  of  these  ftmctions  with  respect  to  v  and  z  are  respectively 
Sc       Sb      /   5y       Sj 

V         -  ... 

"V  (Sv  ~v 


B.  ^  -   C..  ^  =  5  \  J^-^  -  Y^^ 


and 

dC        d3 


V  dz      V  dz      rtz 

If  C  (z)  =  0,  then 

v^    '   -     '  dC        ,   oC 

B   _J1   +B   ^   ;^   =   0, 
V  dz      V  dz  ov 

vhence  —^^0.     Each  term  can  now  be  replaced  by  one  of  the  above  Wronskians,  giving  us 
dz 

rtz    V   V  ov      V  ov  y  dz 

The  quantity  in  parentheses,  hovever,  is  a  constant  multiple  of  the  integral 
appearing  belov;  (  [l]  ^  p.  50);  hence 

_2  .  ^  |i:   r'^K  (2z  sinh  t)  e-2^^dt  =  0. 
nz     n  dz  )o     o^ 

We  shall  restrict  C  to  H  ,  althoi:igh  many  of  the  following  results 

are  valid  for  the  general  cylindrical  function. 

Corollary  6.1.  Let  H  (x)  =0  for  a  and  x  positive.  Then 

a'(x)  >  0  when  a,  >  J  or  x  >  ^;  and  p'(x)  >  0. 

Proof:   Solving  (6.1)  for  -r-  and  separating  real  and  imaginary 
parts,  we  find 

^      =  2x||^|   f  K  (2x  sinh  t)  e'^^"*  cos  2pt  dt 
dx       |dx|  Jo     o 
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2   >uu 

K  (2x  sinh  t)  e'    sin  2pt  dt. 
o  o 


Consider  the  kernel 

f(t)  =  K^(2x  sinh  t)  6"^°-^ 


-2x  sinh  t  cosh  s  -  2at  , 

ds: 
o 


,00 


where 


f(t)  =  J^  e-2^  ^^"^  ^  ^°^^  ^  -  2c^^(-2x  cosh  t  cosh  s  -  2a)ds; 


2 
g(s)  =  (2x  cosh  t  cosh  s  +  2a)  -  2x  sinh  t  cosh  s. 


Since    f(t)  >  0  and  f '(t)  <  0,  the  integral  expression  for  3'  indicates  that 

p'(x)  >  0.   If  we  can  prove  f"(t)  >  0,  it  would  follow  from  Lemma  k.2   that 

a'  >  0.   It  sxiffices  to  consider  g(s). 

2     2  2 

g(s)  =  (2x  cosh  s)   cosh  t  +  (2x  cosh  s)(ikx  cosh  t  -  sinh.  t)  +  ^  ^  ^' 

provided  the  discriminant  of  this  quadratic  in  2x  cosh  s  is  non-positive;  i.e., 

2 

-8a  sinh  t  cosh  t  +  sinh  t  <  0, 

or,  a  fortiori,  a  >  77  .   We  can  show  f  "(t)  >  0  also  when  x  >  |-  ;  for 

2 
g(s)  >  2x  cosh  s  (2x  cosh  t  -  sinh  t) 

=  2x  cosh  s  (2x  sinh  t  -  sinh  t  +  2x)  >  0 

2 
when  1  -  l6x  <  0.  Then  f  satisfies  the  conditions  of  Lemma  ^.2  for  the  above 

values  of  a  or  x. 

Moreover,  with  a  >  0,  we  have  from  (6.I) 
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1  <  2x  |v'  I   Hk^CSx  sinh  t)  dt 

=  4hv|  hJ(x)  h^(x)  =4^  M  IhJ(x)|^ 


where  the  first  equation  results  from  (2.1).     Since 

lEj(x)|<(|;)'^'(l.ij),  C 

as  is  immediately  verified  by  using  Hankel's  asymptotic  formula  (5'5)  with 
Schlafli's  form  for  the  remainder  (3-^))   we  obtain 

Corollary  6.2.   If  H"'"(x)  =  0  with  x  >  0,  then 
l^'WI>   2,V  ,,2  >  V^   • 


B^^ 


Since 


lim   X  |H^(x)|^  =  0  , 
we  note  that 

These  corollaries  have  so  far  been  restricted  to  real  values  x.  For 
complex  values  z,  we  get  the  following  analogue  of  Corollary  6.1,  with  the 
inequalities  suffering  slightly  and  no  longer  independent. 

Corollary  6.3.  Let  r(z)  =  0  with  x,  y,  cl,   P  positive.  Then 

— !-  >  0  when  x  >  7:  and  a  >  "* — ^ —  : 


and 


=<  S;  -  y  *  ^  u  wnen  X  a  2 


^t     *     ^1  >°- 


Proof:  From  (2.2)  for  x  >  0 

(^^     fr.       .   ^.  *\      -2vt,^        f°  r°°-2z  sinh  t  cosh  s  -  2vt 

K     (2z  sinh  t)   e         dt  =  1   e  ds  dt 

^00^  Jo  ^o 


"^    '°°e-P-^^  dt 


Jo    Jo' 
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vhere 

p(t)  =  2x  sinh  t  cosh  s  +  2at 
and 

q(t).  =  2y  sinh  t  cosh  s  +  2pt. 
Sxibstituting  this  in  the  following  form  of  (6.1), 


v'  =  2|zv'|^  nK^(2z  sinh  t)  e"^^*dt, 

ing  real  and  imaginary  parts,  we  get 

a  -  yP  =  2lzv' I   f     e'^  cos  q  dt  ds 
x'^x    '    ^     Jo     )o  ^ 


and 

xP  +  ya  =  2|zv'  I   ]    [   e~^  sin  q  dt  ds 


XX     '    '   ^o  ^o 
where  we  changed  the  order  oT  integration. 

Let  P(t)  =  e~P.   Since  P'  <  0,  q'  >  0  and  q"  >  0, 
we  apply  Lemma  2.2  to  prove  the  second  part  of  the  corollary. 

__      To  prove  the  first  part,  we  shall  apply  Lemma  4.2  to  the  integrand 
P(t)  cos  q(t)  dt,  which  must  first  be  put  into  suitable  form.  We  let  w  =  q(t). 

Since  q'(t)  >  0,  t  =  q"''"(w)  =  t(w) .   Then 


j^PCt)  cos  ,(t)  dt  =  j"  ^^ 


cos  w  dw. 


Consider  the  fiinction 


F(w)    ^ 


:t(v): 


l'  [t(w)]  • 

Then 
and 


q 


•5  Fn(^)  ^  q'(P"q'  -  Pq'«»')  +  5q"  (Pq"  -  P'q')- 
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The  second  term  is  positive.   The  first  term,  apart  from  the  positive  factor 
q'e"^,  satisfies 

g(s)[2y  cosh  t  cosh  s  +  2^]  -  2y  cosh  t  cosh  s  >  [g(s)  -  ij  2y  cosh  t  cosh  s 

when  g(s),  defined  in  the  proof  of  Corollary  6.1,    is  positive,  or  a  >  ^  or  x  >  j-. 

1  2  _/2  -  1 

But  g(s)  -  1  >  0  when  ^  >  -^   and  8ax  +  i+a  >  1  or  a  >  ♦ — ^ —  •   Consequently 

F"(w)  >  0  and  an  application  of  Lenmis  ^.2  completes  the  proof. 

The  inequalities  occurring  in  the  preceding  corollary  can  be  expressed 

more  simply  in  terms  of  the  partial  derivatives  with  respect  to  the  polar 

coordinates  r  and  0  of  the  z-plane. 

Corollary  6.k.      Let  H  (z)  =0  with  x,  y,  a,  p  positive.   Then 

1  Vs  -  1 

If  in  addition  x  >  ^  and  a  >  ^ — ^ ,   then 

Proof:   With  the  help  of  the  Cauchy-Riemann  equations  for  the 
analytic  function  v(z),  we  find 


xP„  +  yP„   xp  +  ya 
3  =  — ^ 2-  =  — >  0, 


r 
and 

o.^  =  -yo-  +  xa  =  -ya  -  xp  <  0. 
P    "^  X     y     X     x 

If  X  and  a  are  restricted  to  the  above  values,  then 

xa  +  ya  xa  -  yP 

r      r  r 

and 

^^   =  -yP  +  xP  =  -yP  +  xa  >  0. 
^0    "^  X     y    "^  X     X 

7.  Behavior  of  the  large  zeros 

In  this  section  we  shall  approximate  the  large  zeros  by  considering 
the  behavior  of  ti  (z)  for  |v|  large  compared  to  1  and  |z|. 


(7.1)     H^(z)  = 


The  first  Hankel  function  may  be  defined  by 
J  (z)    -   e-i^'^J  (z) 


1  sm  v:t 


of  its  limit  if  the  order  is  an  integer.   It  follows  that 


(7.2) 

where 
(7.3) 
and 


,1 


iv  sin  vn   H  (2z)  =  A(v)  +  D(v) 


a/  ^   z       ,  e    ! 

A(v)  =p-(T7y  +-pnT 


D(v)  =  v^ljllJ-- 

k=l 


k  2k 
z 


P  (-v+k+1)  ■*"  p  (v+k+1) 


2  D^(v)  +  D^{v). 


For  V  large^  D(v)  becomes  negligible,  so  we  shall  investigate  the  zeros  of  A(v). 
Then  we  shall  show  that  the  large  zeros  of  E  fcpproximate  those  of  A(v).   Through- 
out, z  will  be  fixed. 

Using  Stirling's  approximation  for  the  gsunma  function,  if  z  =  re  "^ 

iG 
and  V  =  pe   with  0  <  G  <  jt,  we  find  for  large  v  that 

(7-5)     V2^  A{v)^^     B(v) 


exp  (-^  log  V  -   ivn) 


exp  ( V 


If  this  is  to  vanish,  then 


(7.6)     2v  log  —  =  log  1  +  2 


irt 


log  i^  -  "^)  -  exp  (-V  loi 


l^ 


ez  _J 


(2n  +  |)  irt 


for  integral  values  of  n.   Separating  this  into  real  and  imaginary  parts,  we 
obtain 

(7.7)  P  cos  0  log  -2-  -  p(0  -  0)  sin  G  =  0 
and 

(7.8)  2p  sin  0  log  £^  +  (G  -  0)  cos  G  =  (2n  +  |)  jt . 

From  (7.7)  it  follows  that  cos  G  — >  0  as  p  — >  cd  .   Then 


a  =  p  cos  0 


P(|  -  0) 
log|3 


and  making  use  of  the  next  results  we  can  localize  the  large  zeros  as  follows, 
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Theorem  7-1-   Let  H  (z)  =0  vith  0  <  0  <  n.   If  p  is  sufficiently- 
large,  the  zeros  lie  near  the  curve 


a  : 

P(|  -  0) 

3  =  P 

log^ 
^   er 

where  p  satisfies 

2p 

(-  - 
log  P  -  '2 

^^    log 

0f 

P_ 

==  (2n  +  |)  n  . 

As  an  alternative,  we  may  say  that,  for  sufficiently  large  integers  n,  real  and 
imaginary  parts  of  the  zeros  of  H   (z)  are  given  by 

a  =  n(|-   0)    (n.J)    |Tos(n  -:-  1)-|  "^  I  i  +  .     1 

P=.(n.^)    jTog(n.J)|]"'|l.5     I 

where  e     6   =  0  — -i —     . 
n     n        [log  n_| 

In  order  to  show  that  the  large  zeros  of  H  approximate  those  of  B(v), 

we  shall  prove  that  |b(v)|  >  |d(v)|  on  some  circle  centered  at  a  zero  of  B(v) 

and  then  apply  Rouche's  theorem. 

Let  3(v  )  =0,  and  consider  the  circle  C:  v  -  v  =  ee   ,  where 
^  0     '  o 

V  =  pe   and  e  >  0.   Making  use  of  (7-7)  and  (7 .o) ,   \ie   can  show  that 


[/2n  B{v 


)    e"  2 


log  V  +  Ivn 


V   log  v/ez    -V  log  v/ez 
-1  e    >~^     /         +  Q  o  / 


a  -t-  bi    -a  -  bi 
-1  e       +  e 


=  V2  (cosh  2a  +  sin  2b)-'-/^ 


where 


2a  =  2g   cos  CD-  log  -  -  (|  -  0)  sin 
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and 


21)   =  (2n  -  -)   rt+2  €  sin  cd   •    log  —  +  2e  (|  -   0)    cos 
Now  we  get  a  bound  for    |d(v) | . 
-v+2  00 


,    tnK  2k 

(-1;  z 


r    (-V+2)   ^     (k+1)    I    (-V+2) 


euid 


■v+2 


(-v+i)^  (-TT 


IV^)I  = 


-±vn   v+2       00  /    ^  xk  2k 

ve  z  « —  (-1)    z 


r    (v+2)  ,,^     (k+1)    :    (v+2), 


k=0 


-ivn   v+2 
e  z 


"  |(v-i)r  (v) 

for  a  >  0.      Then 

2      /^ 

2  r       A       -V  -ivfl 

As   in  the   case  of  B(v),   we  have 

/Sc  D(v)   e-1/2  log  v+iv« 


<  -^ cosh  a 


k   2r 
Then  to  show  that  |b(v)|  >  |d(v)|  on  C,  it  suffices  to  show  that,  if  c  =  2r  e   , 


cosh  2a  +  sin  2b  >  —^  cosh  a  or 
P 


2        2 
cosh  a  +  sinh  a  >  -sin  2b  where 


-sin  2b  =  cos  2e  sin  o)  •  log  —  +  2€(^  -  0)  cos  co  . 


Choose 


(7-9)     €  = 


2  log 


.^ 


Then  1  -  —  >  -  sin  2b  if  [sin  tn|  >'^  say,  for  p  large  enough.  Now  let  ca   take 
on  the  remaining  values,  i.e.,  | cos  co|  >'^  •     Then 
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(1  -  -2^)  cosh^a  >  (1  -  ^)  cosh^  (|)  >  1 
P  P 

for  p  sufficiently  large.   Consequently  |b(v)|  >  |d(v)|  on  C^  and  Rouche's 

theorem  tells  us  that  the  function  B(v)  +  D(v)  has  as  many  zeros  in  C  as  does 

B(v).  a  similar  argument  would  establish  the  proximity  of  the  zeros  of 

H"""  =  a(v)  +  D(v)  to  those  of  B(v). 

We  note  that  this  proximity  has  already  been  determined  quantitatively 

by  owr   choice  (7-9)  of  e,    if  we  ignore  the  slight  perturbation  introduced  by 

the  approximation  B(v)  of  A(v). 

3.  The  zeros  of  the  Wronskian 

If  the  sphere  is  not  a  perfect  conductor,  the  role  of  H  (ka)  is  taken 
by  the  function 

(8.1)     W^(ka,  k'a)  =  J^(k'a)  ^  H^Cka)  -  H^(ka)  |^  j/k'a) 
where  a  is  the  radius  of  the  sphere  and  k  and  k'  are  the  wave  numbers  of  the 
incident  and  diffracted  waves,  respectively.   By  applying  similar  techniques  to 
W  as  to  H  ,  we  shall  obtain  similar  results.   The  first  of  these  guarantees  the 

V  V 

existence  of  an  infinite  number  of  zeros. 

For  V  large,  since  the  order  of  the  entire  function  m  /    \   is  not 
greater  than  one,  we  have  the  following  relations  for  any  positive  €: 

and 

-r^  J  (z)  =-  J  (z)  -  J  ^(z)   =   0(eP   ). 
dz  v^     z  v^      v+1^ 

1+e 
Moreover,  from  Section  2,  we  know  that  if  y  >  0,  ir(z)  =  oCe*^   )  and  therefore 

n       -,  1+e 

#-H^(z)  =^E^(z)  -  H^^,(z)  =  0(eP   ). 
dz  v^     z     v^  v+1 

Consequently,  the  order  of  growth  of  the  entire  function  W  is  not  greater  than 
lonity  if  Im  k  >  0.  This  is  true  also  if  Re  k  >  0;  however,  the  proof, which 
follows  J  is  more  complicated. 
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Prom  eqxiations  (7-1)  to  (TO),  if  a  <  0,  we  can  easily  obtain 

(ka)     ^  e     (ka. 


„„(..3, 2.-a)~  <'''°'r  (-> 


1    -aP  (-v)  -*"   an(v) 


1+e 


vhence        W     =  0(e^       )    for  a  <  0.     Now  if  a  >  0  and  x  >  0,    then  (see    [l]  ,   p.   27) 
lH^(z)|    =    |c(v)    I(v)| 


where 


and 


1+e 


:(v)    =  0(e^        ) 


|l(v) 


^    cos  ^   "   -/^t  exp[it(-v   +  |)    -   2z   ctg  t] 


.    2v   -)■  1, 
sm  t 


dt 


eyp[pt   -   2x  ctg  t] 


/_  .    2a  +  1, 

o         sm  t 


dt 


<     eP«/2   f       e-2^^(s2  .  -l)-  ds 


<  eP"/2  ^a. 


f  e--   ds    .  f 
•'o  '^l 


-2xs    ,  C~    -2x3     2a, 

ds   +    I      e  s     ds 


where  we  let   ctg  t   ==  s.      The  bracketed  term  is   ].ess   than     yr-  +  ^ ^^  .      Hence, 

(2x) 


H-6 


when  a  >  0  and  x  >  0,  n  (z)  and  its  derivative  equal  0(e    ).   Consequently 

1  +  6 

W  =  0(e    )  for  all  v;    i.e.,  the  order  of  growth  of  V/  is  <  1. 
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Little  now  remains  to  prove  the  next  theorem. 

Theorem  8.1.   If  -  —  <  arg  ka  <  jt  and  -n  <   arg  k'a  <  n, 
then  there  are  an  infinite  number  of  zeros  of  W  (ka.  k'a). 

Proof:  The  function  W^W_^  is  or  order  <  1.  Since  it  ia  an  even 

2 

function  of  v,    it  is  an  entire  function  of  v  ,  with  respect  to  which  its  order 

is  <  — .   Then  it  must  have  an  infinite  number  of  zeros,  provided  we  can  show  it 

is  not  a  polynomial. 

This  may  be  done  by  considering  the  behavior  of  W  W   along  the  p  axis 

by  using  (8.2)  which  follows.   The  same  proof  would  also  indicate  that 
^   2   ,  ,2|   S(rt  -  20) 


W. 


k  -  k' 


for  large  positive  {i,   whence  the  order  of 


W  is  exactly  1,  since  we  know  it  to  be  greater  than  1.   Finally,  we  infer  that 
W.^  /  0  for  lajrge  p. 

To  find  the  distribution  of  the  large  zeros,  we  use  definition  (7-l) 
for  H  .   If  we  let  ka  =  z  and  k'a  =  z'  =  r'e^^,  then  for  |0|  <  it,  |0'|  <  n,  and 
V  not  an  integer,  we  have  for  laa:ge  v 


W  (2z,2z')   = 


la  sm  vrt 


J  (2z') 

V 


z'J'  (2z') 


J   (2z)-e"^^''  J  (2z) 


zJ'  (2z)  -  ze"  ^  J'(2z) 


1 

2ia  sin  vrt 


\l 

\2 

^21 

^22 

where 


and 
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z-2 


21  "  n  (v+i) 

-V 

(-•^-v+i)  ^-^  -  3^^  -  r  (v+i) 


-V  2  -iVTt   V 

A-,0    =     H  /  ..^■^\    (i  -  -tttt)  -  TTTTrrrr    (i  -  tttt);- 


12 


22 


z_ 
v+1' 


z  /        ,    v-2      2v         c  z, 


e  z        /  v+2     2% 

(^    -   TTTT   z    ) 


v+1 


Using  some  elementary  properties  of  the  gcunma  fiJinction  and  the  fact  that 

log  (v+2) '^ log  V  +  — ,  we  obtain 

2v  -ivrt  n  /•      \        'I  o 

(8.2)  -icur  Kh)"  wj2z,  2z')-i  +  ^    "n^.io^^^'^  (-^  -  -'   ) 


r  (v+2) 


1  +  exp 


-  -  2v  log  -  -  log  ^—2 


with  the  aid  of  Stirling's  formula.  Here  we  imposed  the  restriction  0  <  0  <  «  on  v. 
For  the  above  expression  to  vanish,  we  must  have 


2v  log  ^  +  log  g  '^  g  =  (2n  -  -)  irt 


for  integral  values  of  n.  Separating  these  terras  into  their  real  and  imaginary- 
components,  we  find  as  in  the  preceding  section  that  0  approaches  -  as  p  goes 
to  infinity.   Consequently,  we  obtain 

Theorem  8.2.  Let  W  (2z,2z')  =  0  with  0  <  0  <  «,  -rt  <  0  <  it, 
-n  <  0'  <  fl.   If  p  is  sufficiently  large,  then  O-^-^, 

2 

2p  cos  Q  log  ^  -  2p(|  -  0)  +  log  — 2 -'" 

V   ■ 

and 


.2| 


0, 


2p  log  ^   +  2p(|  -  0)  cos  e  -  arg(z^  -  z'^)-^(2n  -  |)  n 
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H^ 


Thus  far,  W  and  H"^  have  had  parallel  treatments  of  their  large  zeros. 
If  ve  consider  large  values  of  z  instead,  by  virtue  of  the  inequalities  obtained 
in  Section  3  'we  know  that  n  (z)  /  0  for  rather  general  conditions.  We  can  prove 
that  W  (z,z')  likewise  does  not  vanish,  but  lacking  any  corresponding  inequalities, 
ve  must  do  this  directly.  We  let  x'  =  z'  =  x'  +  iy'. 

Theorem  8.^.  Let  both  j  z |  and  | z ' |  be  large  compared  to  1  and  | v | ; 


X  >  0,  y  >  0,  x'  >  0,  and  y'  > 


Proof: 


pit 


Then  W  /  0. 


aW 


J.(z-) 


z'J'(z') 


<(z) 


2    ^(^  -  F 

—  e 
jt 


1-1/2    r    I 
z'  '   cos(z' 


rl/2j 


,1/2   .  ,  ,    v:T    7t>     .  1/2 
,'  '   sin(z'  -  p~  -  it)    iz  '  I 


where 


and 


^  /^   .  /    vrt   Its 
■^(z  z'j  '   e  ^    2    4  W 


art 
2 


B  =  y' 


2 


>^     iz  cos  (A+iB)  +  z'  sin  (A+iB) 
=  R  H-  i  I 

>  0 


R  =  sin  A(x  sinh  B  +  x'  cosh  B)  -  cos  A(y  cosh  B  +  y'  sinh  B) 
I  =  sin  A(y  sinh  B  +  y'  cosh  B)  +  cos  A(x  cosh  B  +  x'  sinh  B) . 
Since  the  determinant  of  this  linear  system, 

|(x^  +  x'^  +  y^  +  y'^)  sinh  2B  +  |(xx'  +   yy')  cosh  2B, 


never  vanishes,  neither  does  \1 
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9.   Special  cases 

In  the  preceding  section  ve  let  v  or  z  and  z'  assume  large  values. 
We  shall  now  elininate  this  restriction.   Lacking  any  of  the  integral  representa- 
tions such  as  were  available  for  n,   we  must  content  ourselves  with  much  less 
general  results.  V/e  shall  for  the  most  part  limit  the  variables  to  real  and 
imaginary  values.   This  will  enable  us  to  apply  elementary  considerations, 
particularly  the  fact  that  the  Wronskian  of  linearly  independent  solutions  of 
Bessel's  equation  is  never  zero.   The  first  theorem  is  a  good  illustration  of 
this. 

Theorem  9-1 •   If  z  is  real,  z'  is  positive  or  pure  imaginary,  and 

V  is  real,  then  \-l   (z,z')  /  0. 

Proof:   Using  the  identity  n  (z)  =  J  (z)  +  i  Y  (z),  we  can  rewrite 
the  definition  (O.l)  of  W  as 

d   ,  /,  N      ^  /,  X  d 


W  =  fj  (k'a)  ^  J  (ka)   -  J  (ka)  ~-   J  (k'a)] 
I-  v^     da   V  V     da  v     -' 


+   i  [j  (k'a)  I-  Y  (ka)   -  Y  (ka)  2-  j  (k'a)]  . 
•-  v^     da  v^         v^     da   v^    -' 

Sa.ch  term  in  brackets  is  real,  except  for  a  common  factor  when  k'  is  imaginary 

as  is  evident  from  the  series  representation  for  J  .   If  these  terms  vanish 

simultaneously,  then 

J  (ka)  ^  Y  (ka)   -  Y  (ka)  ^  J  (ke)  =  0. 
v^    da  v^         v^    da  v^ 

But  this  is  impossible.  Recalling  that  z  =  ka  and  z'  =  k'a,  \re  then  obtain  the 
statement  of  the  theorem. 

Theorem  9.21.      If  arg  z  =  —  and  arg  z'  =  0  or  — ,  then  the  zeros 

of  W  (z,z')  occur  in  conjugate  pairs.   If  In  addition  v  =  ip  /  0, 

then  W  (z,z')  ^  0. 

Proof:  Taking  conjugates  of  both  sides  of  (8.l),  when  arg  Ic  =  —  , 
arg  k'  =  0  or  ^,  and  c  =  -e  ^^   or  -1  respectively,  we  find 


This  proves  the  first  part  of  the  theorem.  Ncn^j  when  v  =  ip  ;^  0, 

W   =  ce^''''  [j  (k'a)^H^(ka)   -  H^(ka)  |-  J  (k'a)]  . 

V  l-  -v''     da   v^         v^    da   -v     J 

Then  if  \-I     (and  hence  W  )  vanishes,  the  determinant 

V  V 

j|-J    -J   ^J   =0, 

V  da   -V      -V  da   V 

where  the  argument  is  k'a.   But  this  is  impossible  \-n.th   our  choice  of  v. 

id' 
This  takes  care  of  the  case  in  which  z(=k'a  =  x'  +  iy'  =  r'e  ^  )  is 

real  or  pure  imaginary.   We  wi2J.  now  treat  the  remaining  values  of  z'.   This  will 

permit  us  to  consider,  instead  of  W  ,  a  function  which  is  simpler  in  many  respects. 

The  results  will  be  valid  also  for  real  v. 

Suppose  arg  ka  =  — ;  v  =  ip  or  a;  and  k'  neither  real  nor  imaginary, 

i.e.,  x'y'  /  0.  Then 

W   =  c'  rj-(I?'o)  I-  H^(ka)   -  H^(ka)  ^  J_(k'aO 

V  •-  V      da  V         V     da  v     -^  ' 

where  c'  =  -1  or  -e  ^  according  as  v  =  ip  or  a  respectively.   Conseqnently  W 
(and  hence  W  )  vanishes  only  \ri.th  the  determinant  A(a),  where 

(9.1)     A(a)   =  J  (k'a)  I-  J_(k'a)   -  J-(k'a)  ^  J  (k'a) . 
^  ^        v^     da  v^  V      da   v^ 

Now  we  apply  the   following  forraiila  (Q.J,   p.    20): 

(9-2)  aA(a)      =     (k'^  -  E'^)   jtJ^(k't)   J-(k't)dt, 

whence 

To  show  that  A(a)  never  vanishes,  it  suffices  to  show  that  either  the  derivative 

aA(a) 
or  lya-  — T — -   is  positive,  end  that  the  other  is  non-negative.  V/e  shall  do 

this  on  the  assumption  x'y'  >  0. 

Suppose  first  that  v  =  ip.  Then  from  (9-1)7 

aA(a)  2v  /-k'a/  /Ic'av"^    2   .  ,  ^   -2p0' 

a4?^0+   i-  =  iP  (v^l)r  (-v-^1)  ^"2"^  ^—^      =  -  -  sinh  P«  e 

This  is  positive  if  p  <  0.  Moreover,  the  derivative  in  (9-3)  is  obviously 

non-negative,  since  the  radius  a  is  positive,  aa  always. 
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Now  let  V  =  a.  We  note  that  if  W  =0  then  J  (k'a)  /  0,  for  otherwise 

a  a 

E"''(ka)  =  0,  contradicting  Theorem  5-1.  Then  the  derivative  in  (9-5)  is  positive, 
a 

aA{a) 
Furthermore,  when  a  >  0,  from  (9.I)  we  find  that  the  limit  of  — J— ^  vanishes. 

Then  A(a)  /  0  for  a  >  0,  which  implies  the  same  for  W  . 

A  similar  argument  can  be  made  if  x'y'  <  0,  which  proves 

Theorem  9.3.   Let  arg  z  =  p,  x'y'  /  0,  and  v  =  ip  or  v  =  a. 

Then  W  (z,z')  ?^  0  whenever  Px'y'  <  0  or  a  >  0  respectively. 
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APPENDIX 
Survey  of  resiilts  about  the  zeros  of  H  (x),  x  >  0 

Let  X  be  a  fixed  positive  number.   Let  v  =  a  +  ip  be  a  complex  number 
and  let  H   (x)  be  the  Hankel  function  of  the  first  kind  and  of  order  v  with 

V 

argument  x.  We  have: 

(1)  The  function  H   (x)  is  an  entire  function  of  v  which  vanishes 
for  infinitely  maciy  values  of  v  =  +  v   ,   k  =  1,2,3  *  '  •  v' 

(2)  e'^/^^l^h.)      =  h(1)(x)  ¥  (i-4\ 

k=l  \    V,  / 


where  the  v  depend  on  x. 


(5)  For  values  |v|  which  are  sufficiently  large  compared  to  1  and  to  x, 

the  zeros  of  H   (x)  are  given  by  v  =  +  (a  +  ip  )  where 
V  ^  n   -   n    n 


P,  =  «  (n  .  J)  |lo.  ^n  -  J)  4|  I   E  .  c(ioS^)] 


,(1), 


and  where  n  is  a  (large)  positive  integer. 

This  shows  that  the  distance  of  two  consecutive  zeros  of  H^'^^(x)  tends 

towards  zero,  and  that  the  argument  of  the  zeros  tends  towards  «/2,  although 

their  real  part  tends  towards  infinity.  Therefore,  the  behavior  of  the  zeros 

as  derived  from  the  well-known  formulas  of  van  der  Pol  and  Bremmer  is  not  the 

final  one.   According  to  van  der  Pol  and  Breramer,  we  have  for  large  x  (that  is, 

for  X  »  1),  the  approximate  expression 
1/5 

V  '^  X  +  X  '    T 

n  n 

for  the  n-th  zero  v  of  H^   (x),  where  for  large  n 
n     V 

,2/3 


\~  il^f"  *  11   ^''" 
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and  where 


T   =  1.8Jf2  e 
o 


T   =  5.2^0  e 


i«/5 


;nd  so  on. 


k=l  v. 


(   t^K^(2x  sinh  t)  dt 


(U)   For  the  zeros  v  of  H^  -^(x),  with  k  =  1,2,3  ■•■,   ve  have 

(30 


/ 


K  (2x  sinh  t)  dt 


where  K  denotes  the  modified  Hankel  fimction  of  order  zero. 
o 

(5)   The  following  inequalities  hold:   If  v  =  a  +  iS  and  H^^'^^(av)  =  0, 
where  a  is  real  and  positive  and  ^  >  0,   then  a  <  1. 


x(l), 


and 


and 


If  V  =  a  +  ip  and  H^  ^ {x)    =0  where  x  >  0,  then,  o^  >  0, 
|a|  >  X  -  yx72 


ki   <| 


X  <  P  tan  cut   +  p 


hi  =  i  =4ip|> 


and 


X  ->  0 


a  ->  0,   a  P  -^  0    . 
.(1), 


(6)   If  V  =  a  +  iP  and  H^  ^(x)  =0,  then  a,  P  may  be  considered  as 

functions  of  the  positive  real  variable  x.  V7e  have 

dn,  ^  ^  ,         ^1      ~^  1 
—  >  0  whenever  a  >  tt  or  x  >  t-  . 


dg 

dx 


>  0  whenever  a  >  0  and  x  >  0. 


lim 
X  -)  0 


da   .  dp 
dx     dx 
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